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ABSTRACT 

The average mean wind speed integrated over a disk is shown to be extremely close 
to the mean value of wind speed which would be measured at the center of a disk 
for most geometries in which a WECS (Wind Energy Conversion System) would oper- 
ate. Field test results are presented which compare instantaneous records of 
wind speed integrated over a disk with the wind speed measured at the center of 
the disk. The wind field that a rotating element would experience is presented 
which has been synthesized from the outputs of an array of anemometers. 

INTRODUCTION 

At the present time, the Battel le PNL (Pacific Northwest Laboratories) are en- 
gaged in attempting to statistically characterize the nominal and extreme wind 
fields associated with flow over the disk of rotation of WECS systems. This 
area is presently being addressed experimentally through a field test program 
as well as by theoretical computations. 

In many cases, the only measure of wind speed near a wind turbine may be a sin- 
gle anemometer located at hub height. Questions, with regard to how well wind 
speed measurements from this single anemometer represent the flow field over 
the disk of rotation of the wind turbine and how the measurements relate to the 
wind field a Ijlade element would experience as it rotates, need to be resolved. 

It is the intent of this paper to examine these questions theoretically and 
through the use of field data collected from a set of anemometers arranged to 
describe the swept area of a wind turbine. (See Figure 1.) 

MEAN WIND SPEED 


In general, the mean wind profile in which a wind turbine operates (neutral at- 
mospheric boundary layer) may be described as follows: 

U^ = 1^1nz/Zo 0) 

where = the mean wind speed at height z p = fluid density 

u* = surface shear velocity = A^/p z = height above grade level 

= surface shear stress z^ = surface roughness length 

k = 0.4 
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The ratio of the me^n wind speed averaged of a disk of rotation_of radius R 
and hub height z. , U to the mean wind speed at hub height, U , could be ex- 
pressed as: ^h^ h 



Table 1 gives some example results using Equation (2) to solve for the ratio cf 
the mean velocity averaged over the disk of rotation to the mean velocity mea- 
sured at hub height for various blade radii, hub heights and surface roughness. 

To estimate the ratio of the mean velocity averaged over the disk of rotation 
to the mean wind at hub height, the following formula may be used in lieu of 
Equation (2). 
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The results obtained from Equation (2) indicate that the average wind speed from 
a single anemometer at hub height will be representative of the average wind 
speed over the disk of rotation in almost all cases. However, the instantaneous 
wind speed measured at the hub will not represent the instantaneous average over 
the disk of rotation. Averaging over the disk of rotation will act as a low 
pass filter cutting off much of the high frequency content. The results of a 
field test program performed at Battelle PNL show the effect of special averaging 
over a disk. A schematic of the field test program is shown in Figure 1. 

The field array consists of eight anemometer sets located on the circumference of 
a circle of radius 12.2 m (40 ft) and an anemometer set at the center of the cir- 
cle and at a height of 24.4 m (80 ft). The instantaneous array average (simu- 
lating the average over a disk of rotation) was computed by arithmetically aver- 
aging the instantaneous output from the nine anemometer sets. A 15-second record 
showing both the array average and hub anemometer output is shown in Figure 2. 

The array average is much smoother than the single anemometer trace, as is expec- 
ted. Spectral comparisons of the two traces also show the filtering effect of 
averaging over an area. That is, the spectral amplitude drops off at a quicker 
rate for the disk average than for the single anemometer record. 
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f. 


Using the anemometer array depicted in Figure 1, a wind record (simulating the 
wind field experienced at 12.2 m from the hub by a rotating blade) was synthe- 
sized by incrementally sampling the wind records from the anemometers around the 
circumference of the circle. Each 0.1 second, the wind record was sampled from 
the next anemometer output. 

Thus, each 0.8 second, one complete revolution around the array had occurred. 
This simulates a rotation rate of 75 rpm. A 15-second synthesized wind record 
is shown in Figure 3 along with the wind trace taken from the hub anemometer at 
the same time. The synthesized wind trace shows principally the vertical pro- 
file of the wind (temperature measurements indicated neutral conditions). 

DISCUSSION 

Q. If wind data are available only from a hub-height anemometer, how can the 
full-disk wind flow field be estimated? 

A. In general, the average value of the wind speed over the disk wind flow 
field is extremely close to the average value at the hub, as shown in 
Table 1. However, instantaneously, the two will be different. The aver- 
age oyer the disk could be simulated by filtering the hub anemometer out- 
put with a first-order filter down 3 dB at a frequency of U/D, where U is 
the mean velocity and D is the diameter of the system. Unfortunately, 
filtering processes are inherently time-averaging processes which, if 
used "on line" in a decision loop, have some negative features. 

Q. How do the half-power frequency and turbulence factor change between a 
single anemometer and the total average of all anemometers? How does the 
energy in the 0.1- to 1-Hz range change by using the average of all ane- 
mometers rather than a single anemometer? 

A. In most cases, the special averaging over the disk of rotation does not 
significantly decrease the total amount of energy in the spectrum; thus, 
the half-power point would not change significantly. The changes of the 
half-power frequency would have to be computed on a case-by-case basis. 
Obviously, a small disk does not have as large an averaging effect as a 
large disk. 

The region of the spectra greatly affected by special averaging is the 
high-frequency region. For the atmosphere, this would include the range 
of 0.1 to 1 Hz, the region where most of the energy is lost. Thus, the 
region from 0.1 to 1 Hz would have significantly less energy when area- 
averaged over a disk corrmensurate with a large HECS. Again, the effect 
would be a function of disk size and height above ground level. 
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Battelle welcomes all comments and concerns which identify or assist in identi- 
fying important wind characteristics associated with WECS design. 

1. Formal requests for assistance; Dr. Charles E. Elderkin, Program Manager 
(509) 946-2335 

2. Informal discussion of wind characteristics for WECS design: 

Dr. William C. Cliff (509) 942-5066 

Battelle, PNL 

Department of Atmospheric Sciences 
Richland, WA 99352 


TABLE I. - RATIO OF MEAN VELOCITY AVERAGED OVER A DISK 
OF ROTATION TO MEAN VELOCITY MEASURED 
AT HUB HEIGHT 


Zj^,R Zq (meters) 


(meters) 

0.005 

0.05 

0.5 

5.0 

70:50 

0.99 

0.99 

0.98 

0.97 

50:30 

0.99 

0.99 

0.99 

0.98 

15:10 

0.99 

0.99 

0.98 

0.97 


(a) The terrain features associated with various 
surface roughnesses: 0.005 m = smooth 

0.05 m = moderate 
0.5 m = rough 
5.0 m = extremely rough 
(town) 
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INSTRUMENT ARRAY FOR 
VERTICAL PLANE FIELD PROGRAM 

GILL ANEMOMETERS NO. 1-9 
TEMPERATURE PROBES NO. 10-12 
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Figure 1. - Field test setup for vertical plane test. 


5 


11 


10 


WIND SPEED AT HUB, 24.4m(80ft) 

WIND SPEED AVERAGED OVER DISK 
(RADIUS = 12.2m) (40R) 
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SELECTED DATA ARE STORED IN 9-TRACK 
TAPE AND ARE AVAILABLE TO OTHER TECH- 
NICAL ORGANIZATIONS FOR APPROPRIATE 
APPLICATIONS. TAPES OF DATA SEGMENTS 
MAY BE OBTAINED FROM PNL FOLLOWING A 
DETERMINATION THAT THE DATA ARE 
ACCEPTABLE FOR THE INTENDED APPLICA- 

^'l'S ^ ^ I ^ ^ L 
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Figure 2. - Wind speed at hub versus wind speed averaged over disk. 



Figure 3. - Synthesized wind speed versus wind speed at hub. 




